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I ^ . Abstract. We report on an XMM-Newton observation of the ultraluminous infrared galaxy NGC 6240. The 0.3-10 keV 

spectrum can be successfully modelled with: (i) three collisionally ionized plasma components with temperatures of about 
0.7, 1.4, and 5.5 keV; (ii) a highly absorbed direct power-law component; and (iii) a neutral Fe Ka and K/3 line. We detect a 
' ■ significant neutral column density gradient which is correlated with the temperature of the three plasma components. Combining 

^ ' the XMM-Newton spectral model with the high spatial resolution Chandra image we find that the temperatures and the column 

^ , densities increase towards the center 

' With high significance, the Fe K line complex is resolved into three distinct narrow lines: (i) the neutral Fe Ka line at 6.4 keV; 

I (ii) an ionized line at about 6.7 keV; and (iii) a higher ionized line at 7.0 keV (a blend of the Fe XXVI and the Fe K/3 line). 

fT^ ■ While the neutral Fe K line is most probably due to reflection from optically thick material, the Fe XXV and Fe XXVI emission 

' arises from the highest temperature ionized plasma component. 

, We have compared the plasma parameters of the ultraluminous infrared galaxy NGC 6240 with those found in the local starburst 

CO . galaxy NGC 253. We find a striking similarity in the plasma temperatures and column density gradients, suggesting a similar 

' underlying physical process at work in both galaxies. 

^ : 

^ l' Key words, galaxies: active, AGN - galaxies: individual: NGC 6240 - galaxies: starburst - X-rays: galaxies 

o ■ 

^1. Introduction 1995; Brandt et al. 1997; Kii et al. 1997; Vignati et al. 1999; 

a • Braito et al. 2003, Ptak et al. 2003, Franceschini et al. 2003). 

^ I The Infrared Astronomical Satellite (IRAS) mission 
• '-j . (Neugebauer et al. 1984) resulted in the discovery of a NGC 6240 (IRAS 16504-1-0228) is perhaps one of the most 

large number of galaxies in the local universe (z < 0.3) with interesting and best studied objects in this class. Although 

extraordinarily large infrared luminosities (Lir > 10'^ Lq for it is slightly less luminous (Lfir - 6 ■ lO" Lq for Hq - 

■ ■ ' Ho = 75kms"^Mpc '; see Sanders & Mirabel 1996 for a 75 km s ' Mpc ' ; Genzel et al. 1998) than the typical IR lumi- 
review). The number density of these Ultraluminous InfraRed nosity used to define a ULIRG, it exhibits all of the ULIRG 
Galaxies (ULIRGs) exceeds that of optically selected Seyfert characteristics, and it is thus studied as a prototype of this 
galaxies and QSOs with comparable bolometric luminosities class of sources. Optical/IR observations show a complex sys- 
(Soifer et al. 1987; Sanders et al. 1988a & 1988b; Sanders & tem with a disturbed morphology consisting of two gravitation- 
Mirabel 1996) by a factor of ~ 1.5-2 (Sanders et al. 1999). ally interacting nuclei (Fosbury & Wall 1979; Fried & Schulz 
Some of the most fundamental issues regarding ULIRGs 1983). An 750 observation of NGC 6240 (Genzel et al. 1998) 
are: the nature of the power source (strong starburst and/or a favors a starburst-dominated power source, in which the high 
dust-enshrouded AGN); and the relative contribution to their infrared emission arises from warm, absorbing dust, surround- 
bolometric luminosity. The capability of X-rays to penetrate ing the inner parts of the galaxy. Early X-ray observations 
high absorbing material makes X-ray observations essential of NGC 6240 with ROSAT (Schulz et al. 1998) revealed an 
in attempting to understand the physical processes at work in exceptionally high X-ray luminosity (Lo.i_2.()keV ~ IO'^Lq), 
ULIRGs. A considerable portion of ULIRGs have been found which further confirmed the importance of NGC 6240 in un- 
to contain a hard X-ray source, highly absorbed by a molecular derstanding the AGN-starburst connection. Therefore, consid- 
torus, which indicates the presence of a hidden AGN (Mitsuda erable effort has been made to analyze the spectral and spa- 
tial properties using different X-ray satellites. Moreover, NGC 
6240 can be considered as the prototype objects of highly db- 
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First ROSAT observations showed that a hot, diffuse gas 
model in combination with a power-law component can give 
a good description of the (0.1 - 2.4) keV energy band spectrum 
(Schulz et al. 1998). The power-law component was interpreted 
as evidence for a hidden AGN in NGC 6240. A luminous ex- 
tended soft X-ray component was reported by Komossa, Schulz 
& Greiner (1998) using ROSAT H^l observations. This emis- 
sion should be responsible for at least 60% of the 0. 1 - 2.4 keV 
luminosity (within a region of 5" radius from the center). The 
combination of the spatial extent and the Raymond-Smith-like 
thermal spectrum led to the idea that the origin of the soft X- 
rays were strong thermal processes indicative of a powerful, 
X-ray luminous starburst. 

According to data from ROSAT and ASCA, Iwasawa & 
Comastri (1998) reported the detection of a highly absorbed 
hard X-ray source in NGC 6240. They pointed out that the soft 
X-ray spectrum (~ 0.5 - 2.0 keV) could be described by a su- 
perposition of two optically thin thermal plasmas with different 
temperatures (kTi ^ 0.2-0.6 keV and kTa - 1 keV) and an Nh 
value roughly an order of magnitude higher than the Galactic 
value of Ne.gai - 5.8 ■ 10^** atoms cm"^. At higher energies 
(~ 3 - lOkeV) ASCA observed an extremely hard continuum 
slope accompanied by the presence of a strong Fe Ka com- 
plex (Mitsuda 1995; Iwasawa & Comastri 1998). This can be 
interpreted in terms of a luminous AGN completely hidden by 
Compton-thick material, whose emitted light is reflected into 
the line-of-sight (e.g. George, Nandra & Fabian 1990). In ad- 
dition, the Fe complex could be well explained by a superpo- 
sition of two Gaussian lines at Ei = (6.44 + 0.04) keV and 
E2 - (6.87 ± 0.05) keV, this latter was interpreted as a strong 
hint for a highly ionized reflecting material. 

Vignati et al. (1999) and Ikebe et al. (2001) investigated 
NGC 6240 by using data from BeppoSAX and RXTE, respec- 
tively. Both authors assert the existence of an AGN component 
directly related to a Compton thick absorber (with Nh - 2 • 10^"* 
atoms cm"^). 

Recent Chandra observations have spatially resolved the 
emission from the two nuclei and the surrounding starbursts at 
X-rays. The presence of the neutral Fe Ka lines and the hard 
power-law-like X-ray emission from both nuclei was used to 
infer that both nuclei of the galaxy harbours an AGN (Komossa 
et al. 2003). In addition, they found residua in their spectral fits 
at energies above the neutral Fe K line, which they interpreted 
as possible emission from H-like Fe from the nuclear regions. 

This paper contains the analysis of the X-ray data of 
NGC 6240 taken with XMM-Newton (Jansen 1999). The data 
analysis and details of the XMM-Newton observations are de- 
scribed in Section 2. Section 3 contains the spectral fitting re- 
sults. A comparison with the Chandra observations is given 
in Section 4. The summary of the XMM-Newton results on 
NGC 6240 is given in Section 5. We use a Hubble constant 
of Ho - 50 km s"^ Mpc"^ and qo = throughout. 

2. Technical details: XMM-Newton observation and 
data processing 

NGC 6240 was observed with XMM-Newton for the first time 
during orbit 0144 on the 22"'' of September, 2000, with aU in- 



struments functioning except for the RGS. The EPIC pn camera 
was operated in the full-frame mode and the EPIC MOS cam- 
eras were operated in the large-window mode. Each of the de- 
tectors used the medium filter The exposure times were 25150 
s and 24349 s for the pn and MOS, respectively. The second 
observation, with a shorter exposure time of about 18 ks, was 
performed during orbit 0413 (12* of March, 2002), with afl in- 
struments working in the same mode and with the same filters 
as in the first exposure. In addition, the RGS was functional and 
operated in the Spectro-i-Q mode, however the exposure time 
was too low to collect sufficient signal (raw exposure times: 
18917 s for each of the RGS devices). 

The MOSl and MOS2 data sets from each observation 
were merged for this analysis. Single and double events have 
been selected and the resulting PHA file was grouped with 
a minimum of 20 counts per bin. After correcting the MOS 
data for background flaring events, the GTI corrected exposure 
times for MOS 1 and MOS2 were 2338 1 s and 23601 s (first ob- 
servation - orbit 1 44) and 11318s and 1 1 60 1 s (second obser- 
vation - orbit 0413). This gives us a total of about 35 ks worth 
of exposure time which is used in the analysis. Examining each 
observation separately and comparing them to the merged data 
set, we found no discrepancies in the spectral fitting parameters 
within the data uncertainties, justifying the merging process. In 
addition, the EPIC pn data sets were merged in a similar fash- 
ion (orbit 0144: 16255 s, orbit 0413: 9750 s). 

The EPIC MOS energy resolution (FWHM) is slightly bet- 
ter than that of the EPIC pn. Comparing Fig. 25 and 26 of 
the XMM-Newton User's Handbook (v2. 1 ; Document number: 
XMM-PS-GM 14') one gets an energy resolution for the EPIC 
MOS cameras of about 140 eV between 6 and 7 keV. For the 
EPIC pn camera the energy resolution ranges between 160 and 
240 eV, but it is better for double events (upper curve of Fig. 
26) and decreases for single events when considering on-axis 
sources. As we have used single and double events, the EPIC 
pn energy resolution ranges between 160 and 240 eV (the exact 
value depends on the ratio between single and double events 
which is presently not available). In order to utilize as many 
photons as possible we have concentrated on the EPIC MOS 
data for subsequent spectral analyses (about twice as many 
photons than in the EPIC pn data). With the combined MOS 1/2 
photons from the two observations we collect in the 0.3-10 
keV energy range about 1.0 ■ 10"^ photons cm"^ s"' keV~' 
(compared to 5 • 10 ''photons cm"^ s ' keV ' for EPIC pn). 
Nevertheless we use the EPIC pn data to independently check 
the EPIC MOS results for consistency. 

The data were processed using the standard procedures 
of the XMM-Newton Science Analysis System (XMM-SAS^) 
version 20020413^031-5.3.0 with default parameters for the 
EPIC chains. These event lists were calibrated with the latest 
available calibration files^. 

Source and background counts were extracted from circu- 
lar regions with radius of 70". The background photons were 
extracted from a source free region with the same radius. For 

' ftp : //xitiiti . vi 1 spa . esa . es/p ub/A03 /XMH_UHB . pdf] 
^ http://xiiim.vilspa.esa.es/sas/ 
^ .http : //ximti ■ vilspa . esa . es/ccf /^ 
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Fig. 1. The (0.3 - 10) keV spectrum of the combined source and 
background counts (upper curve) and the background counts 
(lower curve with the black dots) for the EPIC MOSl/2 data. 
The spectrum is not background dominated below 10 keV and 
the Fe K line complex can be investigated with a high signal- 
to-noise ratio. 

the X-ray spectral analysis we have used XSPEC version 1 1.2 
(Arnaud 1996) as well as FTOOLS version 5.2. The quoted er- 
rors on the derived best-fitting model parameters correspond to 
a 90% confidence level. 

3. Spectral analysis 

The combined MOSl/2 spectrum of NGC 6240 is not back- 
ground dominated below 10 keV (c.f. Fig. 1); therefore, the en- 
ergy range between (0.3 - 10)keV is analyzed with relatively 
good signal-to-noise. 

3.1. The soft energy spectrum 
3.1 .1 . Line identification 

First, we searched for emission lines in the soft spectral re- 
gion (0.3 - 3.5 keV). Figure 2 shows the EPIC MOSl/2 source 
spectrum and several identified emission lines. For illustration 
purposes we have included a 1-keV-model spectrum of a dif- 
fuse collisionally-ionized plasma, calculated using the APEC 
algorithm (Astrophysical Plasma Emission Code"^; Smith et al. 
2001) with variable abundances as provided by the vAPEC 
model in XSPEC (solid line). The abundances of O, Ne, Mg, 
Si, S, and Ar were fixed to 10 times solar and the abundance 
of Fe was fixed to 0.01 times solar (the others were fixed to so- 
lar abundances). Pairs of emission lines, which are in different 
ionization states, can be resolved in the EPIC MOSl/2 spec- 
trum. Moreover, higher energy emission line such as S XV and 
Ar XVll are also detected. The ratio of the emission line pairs 
and the emission lines at higher energies enable us to deter- 
mine how many thermal plasma components are required to fit 
the soft energy spectrum. 

.http : //hea-www . harvard . edu/APEC/ 1 
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Fig. 2. EPIC MOSl/2 spectrum of NGC 6240 (upper 
curve). For illustration purposes a vAPEC model with a 
temperature of 1 keV has been included (lower line). 
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Fig. 3. Fitting only one vAPEC model to the EPIC MOSl/2 
data in the (0.3-2.5) keV range, where O, Ne, Mg, Si, S; the 
other a elements; and Fe-Ni abundances have been left free. 
The resulting fit is statistically unacceptable. 



3.1 .2. Tine presence of two distinct plasmas 

Applying only one vAPEC model to the EPIC MOSl/2 data 
with thawed abundance parameters for O, Ne, Mg, Si, S, the 
a elements, and Fe-Ni, did not result in an acceptable fit. As 
can be seen in Fig. 3, significant residua remain in the Fe L 
emission line region around 0.8 keV, as well as in the (1.0 - 
1.5) keV range (x^ = 1.5). The observed Fe L bump, the flux 
ratio between two emission lines in different ionization states 
and the fit residuals all suggest the presence of a second thermal 
component. 

We have added a second vAPEC component to the first 
model, which resulted in an acceptable fit (x^ = 10). The sig- 
nificance of this new component was tested with the F-test. The 
second component is required with a significance greater than 
5cr. The pn data explicitly confirm the EPIC MOSl/2 results. 
The best fit parameters of the two vAPEC components are re- 
ported in Table 1 . 
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3.2. The Fe K line complex 

The high throughput of XMM-Newton allows for precise mea- 
surements of the Fe K line complex compared to previous X- 
ray missions. In this Section we will demonstrate with high 
significance that the Fe K line complex is resolved into three 
distinct Unes. The line energies are given in the rest frame of 
the source. 

3.2.1 . XMM-Newton detection of three distinct Fe K 
emission lines 

As a simple first approach the 5 -7.2 keV energy range has been 
fitted with a power-law model and two Gaussian lines (c.f. Fig. 
4, upper panel). The centroid line energies are found at (6.5 
± 0.3) keV and (6.68 + 0.02) keV. Strong residua remain at 
around 6.9 keV and the fit is statistically unacceptable (x^ = 30 
for 22 d.o.f.). We can exclude the possibility that the Fe line 
complex can be modelled by a narrow unresolved 6.4 keV fine 
and a blend of higher ionized lines (Fig. 4, middle panel; = 
30 for 20 d.o.f.). 

The strong residuals which remain at 6.8 - 6.9 keV clearly 
are reduced by using a power-law plus three narrow unresolved 
Gaussian lines (x^ = 15 for 19 d.o.f., see Fig. 4, bottom panel). 
The centroid line energies and corresponding 90 per cent co- 
incidence ranges are (6.41 + 0.02) keV, (6.68 ± 0.02) keV and 
(7.01 + 0.04) keV. With high statistical significance, the XMM- 
Newton data reveal for the first time the presence of three nar- 
row Fe K lines. 

3.2.2. Tlie possible physical origin of the ionized Fe K 
lines 

One plausible model to explain the presence of ionized Fe K 
emission lines arises from the presence of a collisionally ion- 
ized plasma (using the APEC model in XSPEC). The plasma 
with a temperature of about 5.5 keV produces the Fe XXV (at 
6.68 keV) and Fe XXVI (6.93 keV) emission fines (c.f. also 
the discussion of the broad band spectrum in Section 3.3). In 
the spectral fit for the Fe line complex we have also included 
the emission from the Fe Ka line at 6.4 keV and the Fe K/3 fine 
at 7.058 keV, each modelled with a Gaussian fine. The normal- 
ization in the spectral fit of the line intensity of the Fe K/3 line 
was set to 1/8.8 of the Fe Ka fine intensity, as expected from 
atomic physics. Fig. 5 shows the EPIC MOSl/2 spectral fitting 
results for the Fe line complex in NGC 6240 (Xr = 1.1). The 
model components are displayed in Fig. 6. 

In addition, we tested other possibilities to describe the ori- 
gin of the ionized line emission in NGC 6240. First, we mod- 
elled the Fe K line complex by He- and H-Uke emission from 
the accretion disk. Using the Ballantyne et al. (2001) model 
available within XSPEC we are unable to find a statistically ac- 
ceptable fit. This is expected as BeppoSAX observations show 
that the direct AGN component is heavily absorbed and only 
dominates the spectrum above 8 keV (c.f. Fig. 3 of Vignati et 
al. 1999). Therefore, the putative He- and H-Uke emission from 
the accretion disk should remain undetectable. 



Secondly, we substituted the APEC model by an absorbed 
power-law component plus three Gaussian lines. The fit to the 
(0.3 - 10)keV spectrum was statistically unacceptable (x} - 
1.8) and is unable to explain the broad-band spectral energy 
distribution, nor the presence of ionized Fe K Unes. 

Thirdly, fitting the broad-band spectrum with only two 
APEC components was statistically unacceptable. A third 
APEC component is required by the data. The two lower tem- 
perature APEC components have already been discussed in 
section 3.1.2. 

Colbert et al. (2002) have modelled the Fe K features in 
NGC 1068 from both optically thick and optically thin ionized 
reflection emission. The optically thick case was modelled us- 
ing the PEXRIV model in XSPEC, the optically thin case by a 
power-law plus edge model. The authors find that the Compton 
reflection component provided by the optically thick model 
does give a better fit than the optically-thin reflection model. 
We have modelled the XMM-Newton spectrum by substituting 
the APEC model wifli the PEXRIV model (and by adding two 
Gaussian lines to model the ionized Fe K Unes). This model 
does not provide an acceptable fit (x^ = 4.6). The same holds 
for the power-law plus edge model for an optically thin warm 
reflector. The XMM-Newton data do not require, within the 
available statistics, an additional ionized edge. 

The collisionally ionized plasma provides the best spectral 
fitting results and appears the most likely explanation for these 
XMM-Newton observations (see Table 1 for the best fitting pa- 
rameters). 

3.3. The broad-band X-ray spectrum 

The best fit to the 0.3 - lOkeV XMM-Newton spectrum of 
NGC 6240 is obtained using: (i) three dififerent plasma tem- 
peratures; (ii) a highly absorbed power-law component; and 
(iii) a neutral Fe Ka and K/3 line (c.f. Table 1 and Fig. 7). 
The broad-band fit with these spectral parameters results in an 
acceptable fit (x^ = 0.93 for 190 d.o.f.). In the soft energy 
band the spectrum is dominated by the cold (kT = 0.7 keV) 
and medium (kT = 1 .4 ke V) plasma components. The medium 
component is intrinsically absorbed with an absorbing column 
of (3.6 + 0.3) • 10^' cm"^. For the a elements we find el- 
ement abundances of (0.15 + 0.08) relative to solar values. 
The Fe and nickel abundance is lower compared to that of 
the a elements with (0.05 + 0.01); however, not significantly 
different. The cooler plasma component has an intrinsic col- 
umn density of (2.0 + 0.3) • 10^' cm"^, and a temperature 
of kT = (0.66 ± 0.03) keV. The available photon statistics 
does not allow us to constrain the element abundances for the 
cooler plasma component; therefore, we have fixed the a and 
Fe and nickel abundances to that of the hotter plasma compo- 
nent, which results in a statistically acceptable fit for the soft 
energy region. However, this should be handled with caution. 
Finally, including a third hot plasma component and an ab- 
sorbed power-law component allows us to explain the origin of 
the ionized, Fe XXV and Fe XXVI line emission which were 
discussed in detail in Section 3.2.2. This hottest plasma compo- 
nent has a temperature of kT = (5 .5 ± 1 .5) keV and an intrinsic 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




vAPEC(l) 


vAPEC(2) 


APEC(3) 


power-law 


FeKo- 


FeXXV 


Fe XXVI 


FeK/3 


NH(fit) [10-2 cm-2] 


0.20 ± 0.03 


0.4 ± 0.03 


4.1 ± 1.3 


100 ± 20 


_ 


_ 


_ 


_ 


kT [keV] 


0.66 ± 0.03 


1.4 ±0.2 


5.5 ± 1.5 


_ 


_ 


_ 


_ 


_ 


Norm " 


4.8 ± 1.2 


2.5 ± 1.5 


3.7 ± 1.3 


6.0 ± 3.0 


2.6 ± 0.8 


1.8 ±0.7 


0.9 ± 0.7 


0.3 ±0.1 


f [10-'2ergcm-2 s-'] 


1.5 ±0.4 


0.7 ± 0.4 


1.4 ±0.5 


40 ±20 


0.25 ± 0.08 


0.19 ±0.07 


0.10 ±0.08 


0.03 ±0.01 


L [10^*2 ergs-'] 


3.9 ± 1.0 


1.9± 1.1 


3.6 ± 1.3 


100 ± 50 


0.60 ± 0.20 


0.50 ± 0.20 


0.30 ± 0.20 


0.08 ± 0.02 


Equivalent width [eV] 










300 ± 100 


220 ± 90 


120 ±90 


50 ±20 



Table 1. Spectral fit components to the (0.3-10) keV energy band. The individual columns 2 to 8 refer to the spectral compo- 
nents applied in the fit. The rows specify the measured values. The broad-band spectral energy distribution measured with these 
components is shown in Fig. 7. The foreground absorption has been fixed to the Galactic value. ° NOTE: For the first three 
components the normalization is given in units of [10-'^/(47r(DA ■ (1 + z))^) Jng ■ ■ dV], where Da is the angular size distance 
to the source (cm), ng is the electron density (cm^^), and ne is the hydrogen density (cm^^). For the last five components the 
normalization is given in units of [10^^ photons cm"^ s ' keV '] at 1 keV. The fluxes and luminosities are absorption coiTected. 



absorption of (4.1 ± 1.3) ■ 10^^ cm ^. The APEC components in 
this model are regarded as a description for starburst activity. 

A heavily absorbed power-law is required by the XMM- 
Newton data. If neglected, significant residua remain above 
about 7 keV and the fit is not statistically acceptable. We can 
confirm the BeppoSAX power-law with the XMM-Newton data 
at a 4 cr significance level (according to an F-test). The absorb- 
ing column density of (1.0 + 0.3) • lO^'* cm~^ is a factor of two 
smaller than that found by BeppoSAX. We believe that the pho- 
ton index can be better measured by BeppoSAX gwtn the larger 
energy range over which it was determined. Therefore, we have 
fixed the power-law photon index to 1.8. The absorbed power- 
law is indicative of the direct emission from the accretion disk 
corona. 

We also note that with the present statistics the fit is unable 
to independently determine the temperature of the hot plasma 
component, and the power-law photon index, and the column 
density. 

In order to cross-check the EPIC MOS results obtained 
from the (v)APEC modelling with a different plasma code 
we substituted the three components with the coiTesponding 
MEKAL model: in the soft regime the two vMEKAL models 
and in the hard part of the spectrum a MEKAL model. We find 
no significant difference in the absorbing column densities and 
plasma temperatures with respect to the values listed in Table 
1 . For the Nh values we get (0. 15 ±0.02; 0.5 +0.3; and 3 .5 + 1 . 1 ) 
•10^^ cm"^. The corresponding plasma temperatures are kT = 
(0.62 + 0.03) keV, (1.2 + 0.4) keV and (6.0 + 2.6) keV. 

A second test was performed to compare the results of the 
combined EPIC MOS data with the data of EPIC pn. The re- 
sulting values Nh - (0.19 + 0.03) ■ 10^^ atoms cm"^ with 
kT = (0.58 + 0.03) keV, and Nh = (0.4 + 0.1) ■ lO^^ ^^^^^ 
cm-2 with kT = (1.2 + 0.3) keV for the two vMEKAL com- 
ponents are in good agreement to the values of the EPIC MOS 
data. The fitting with a MEKAL model yields in the values 
Nh = (3.8 + 2.0) ■ lO^^ ^toms cm^^ and kT = (5.0 + 1.9) keV. 

3.4. Comparison with local starbursts 

A striking similarity in the plasma parameters and absorbing 
columns is found between NGC 6240 and the prototype lo- 



cal starburst galaxy NGC 253. Pietsch et al. (2001) have an- 
alyzed the nuclear emission of NGC 253 with XMM-Newton. 
They found that the nuclear spectrum can be modelled by a 
three temperature plasma with the higher temperatures increas- 
ingly absorbed. The derived temperatures for the nuclear re- 
gion are 0.56, 0.92 and 6.3 keV. The coiTesponding absorbing 
column densities are (0.34, 1.78, 13.2)-10^^ atoms cm"^. These 
temperatures and the corresponding absorbing column densi- 
ties are remarkably similar to those found in NGC 6240. In 
addition, Pietsch et al. detected He-like Fe K emission at 6.68 
keV, similar to the detection in NGC 6240. The non-detection 
of the weaker H-like line in NGC 253 might be due to a lower 
signal-to-noise compared to NGC 6240. The authors attribute 
the presence of the high temperature component to emission 
from type lb and type Ila SNRs. In addition, similar soft X-ray 
emission lines from the two cooler plasma components are de- 
tected in NGC 253, e.g. emission from Mg XI, Mg XII, Si XIII, 
Si XIV (c.f. Section 4. 1 of Pietsch et al). The presence of differ- 
ent plasma components in the prototype local starburst galaxy 
NGC 253 support our spectral modelling, and further suggests 
that there may be a common underlying physical mechanism, 
which appear to be dominated by emission from type lb and 
type Ila SNRs. 

4. Comparison with Chandra 

In the following we combine the high throughput of XMM- 
Newton with the high spatial resolution of Chandra. While 
Chandra has resolved the nuclear region into two nuclei and 
diffuse emission components (Komossa et al. 2003), XMM- 
Newton allows for X-ray spectroscopy with high statistics. 

The inner 30" remain unresolved with XMM-Newton as ob- 
tained from the PSF The PSF was created simulating a point 
source at the same position on the detector as NGC 6240, and 
at an energy centered at 6.5 keV, using the XMM science simu- 
lator software SciSim^. The energy selected coiTesponds to the 
mean energy of the Fe K line complex. Serendipitous sources 
(locaUzed through the XMMSAS task emldetect) and out of 
time events were removed from the data. Each photon in the 
profile was corrected for vignetting and divided by the expo- 

^ |http : //ximti ■ vilspa . esa . es/scisim/, 
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sure map, to remove chip gaps and hot pixels. Since neither the 
simulated source nor NGC 6240 show any pile-up, the sources 
were normalized to the maximum flux. Using such a PSF will 
allow us to further investigate the spectral components around 
the two nuclei detected with Chandra. The data were accumu- 
lated in the (5.5 7.5)keV energy range for the combined 

MOSl/2 data. The radial intensity profile of the emission, to- 
gether with the PSF are shown in Figure |8] Due to the large 
width of the PSF, the (PSF-uncorrected) extended emission ap- 
pears to start at a distance greater than 30" from the center. The 
spectral analysis of the extended emission was extended up to 
a distance of 70". Beyond this distance the number of counts is 
too small to make a significant contribution to the statistics. 

We have performed a spectral analysis of the inner 30" in- 
dependent from the results obtained from the total emission 
(0-70"). Interestingly, the spectral fit to the inner 30" is consis- 
tent with the spectral fit components and parameters obtained 
for the total emission (c.f. Fig. 7) except for the flux emitted 
by the cold APEC component. Ninety-eight per cent of the 
(0.3 - 10 keV) flux is emitted in the inner 30" core. The derived 
temperatures for the plasma components confirm the existence 
of different plasma contributions in the core region. 

In the following we compare the XMM-Newton plasma 
components with the energy bands used in the high resolu- 
tion Chandra image (Fig. 3 of Komossa et al. 2003). In the 
softest energy band (red coloured; 0.5 -1.5 ke V) the emission 
from the 0.7 keV plasma component dominates the spectral en- 
ergy distribution, with a contamination level of about 10 per 
cent from the 1.4 keV component. In the (1.5 - 5)keV band 
(yellow and green coloured), we have a mixture of all three 
plasma components and cannot uniquely attribute a single tem- 
perature to the spatial location. Above 5 keV the hot 5.5 keV 
plasma component dominates the spectrum and is responsible 
for the emission from Fe XXV and Fe XXVI. This is confirmed 
by the Chandra observation; Komossa et al. (2003) note that 
the neutral Fe K line dominates the spectral energy distribu- 
tion between (6-7) keV, and noticed that some residua remain 
in their spectra which can be attributed to ionized Fe K emis- 
sion. In addition, the XMM-Newton data reveal a strong column 
density gradient for the diff'erent plasma components for the 
inner region, decreasing outwards from about 4.1 ■ 10^^ cm"^ 
(kT = 5.5 keV), to 0.4 ■ 10^^ cm^^ (kT = 1 .4 keV), and down to 
0.2 ■ 10^2 cm-2 (kT = 0.7 keV). 



5. Summary 

XMM-Newton observations of the ULIRG NGC 6240 reveal: 
(1) For the first time, the presence of three distinct Fe K lines 
located at (6.41 + 0.02) keV, (6.68 + 0.02) keV and (7.01 + 
0.04) keV (in the rest frame of the source); 

(2) That the broad-band spectrum is successfully described 
by the emission from three diff'erent plasma components with 
temperatures of (0.66 + 0.03) keV, (1.4 + 0.2) keV, and (5.5 
+ 1.5) keV. The ionized lines can be explained by emission 
from the highest temperature plasma components resulting in 
emission from Fe XXV and Fe XXVI. The three plasma com- 
ponents are found at diff'erent locations in the Chandra image. 



showing that the temperature and the column density are in- 
creasing towards the center of NGC 6240; 

(3) A strong column density gradient for the different 
plasma components increasing with temperature from (0.20 + 
0.03) -1022 cm-2, to (0.40 + 0.03) -lO^^cm^^ to (4.1 + 1.3) 
■1022 cm-2; 

(4) The presence of several emission lines in the soft en- 
ergy band which is in agreement with previous ASCA mea- 
surements; 

(5) A highly absorbed direct power-law component, con- 
firming the BeppoSAX results; 

(6) A striking similarity to the local starburst galaxy NGC 
253. In both, NGC 253 and NGC 6240, three plasma compo- 
nents are required, and their temperatures and column density 
gradients are remarkably similar; hence suggesting a similar 
underlying physical process at work in both galaxies. 

The spectral capabilities of XMM-Newton in combination 
with the high spatial resolution of the Chandra telescope gives 
us a more precise understanding of the nuclear region in NGC 
6240 (c.f. Fig. 9). 
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Fig. 9. We have binned the Chandra data in three different energy bands as seen by the XMM-Newton spectrum. The crosses 
mark the position of the two nuclei as detected with Chandra. The left panel shows emission from the neutral Fe K line (selected 
between 6.1 and 6.35 keV in the observers frame, emission from the ionized lines (middle panel, energy range between 6.4 and 
7 keV) as well as emission between 7 and 10 keV from the highly absorbed power-law component (right panel). The spectral 
capabilities of XMM-Newton in combination with the high spatial resolution of the Chandra telescope gives us a more precise 
understanding of the nuclear region in NGC 6240. 
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Fig. 4. Upper panel: Spectral fit to the XMM-Newton com- 
bined MOS 1/2 data of NGC 6240 in the Fe K Une range. Fitting 
the Fe line complex with only two narrow lines and a power- 
law continuum results in a statistically unacceptable fit. 
Middle panel: Fitting the Fe line complex with a narrow neu- 
tral Une and one broad line also result in a statistically unac- 
ceptable fit. 

Bottom panel: An acceptable fit is obtained when fitting three 
Gaussian lines above the power-law continuum. The centroid 
line energies are at (6.41 ± 0.02) keV, (6.68 ± 0.02) keV and 
(7.01 ± 0.04) keV. With high statistical significance, the XMM- 
Newton data reveal for the first time the presence of three nar- 
row Fe K Unes. 




5.5 6 6.5 7 

channel energy (keV) 

Fig. 5. Spectral fit to the XMM-Newton MOS 1/2 data of NGC 
6240 in the Fe K line range (with the identifications of the Fe 
lines). Emission from a collisionally ionized plasma has been 
assumed to model the ionized Fe K line emission. The neutral 
6.4 Fe Kor line has been modelled by a simple Gaussian line. 
We have also included the emission from the Fe K/3 line where 
the ratio has been fixed to 1/8.8 according to atomic physics. 
The Fe XXV and Fe XXVI Unes originate from the hottest 
plasma component. 



n 



J 



- 6 6.2 6.4 6.6 6.8 7 7.2 

Energy 

Fig. 6. Spectral model for the Fe K line complex. Two unre- 
solved Fe lines (Ka and K/5) and emission from a collision- 
ally ionized plasma (APEC model in XSPEC) can explain the 
emission from the Fe K Une complex in NGC 6240 (for Une 
identifications see Fig. 5). 
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channel energy (keV) 

Fig. 7. Spectral fit to the merged MOS 1/2 data of NGC 6240 in 
the 0.3-10 keV energy range. In the soft energy band, a statis- 
tically acceptable fit is obtained by using a hot and cool plasma 
component, where the a and Fe-Ni abundances have been left 
free in the fit. The emission from a third hot plasma compo- 
nent is also required to reproduce the line ratio of the ionized 
line emission of the Fe K complex. A narrow unresolved Fe Ka 
and Fe K/3 line, with a relative ratio of 8.8 : 1, as expected from 
atomic physics, has been added to model. A highly absorbed 
(Nh = (1.0 + 0.3) ■ lO^'* cm"^) power-law component with an 
photon index of 1.8 significantly improves the fit above about 
7 keV. The fit to the 0.3-10 keV band is statistically accept- 
able with the above mentioned components (x^ = 0.93 for 190 
d.o.f). 




Fig. 8. The radial intensity profile of the EPIC MOS 1/2 data 
compared to the instrumental point spread function (PSF; nor- 
malized to the maximum flux). The solid line shows the MOS 
PSF of 6.5 keV and the photons have been selected in the 5.5- 
7.5 keV range to match the emission from the Fe line complex. 
Extended emission is detected at distances of about 30" from 
the X-ray centroid position. 



